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Abstract. Using new techniques based on the polarimetric fluctuation spectrum, the fluctuation behaviour of the polarization 
of individual pulses is examined in three pulsars that show drifting subpulses, allowing various aspects of the fluctuations to 
be quantified for the first time. Of the three pulsars studied, only PSR B0809+74 shows behaviour completely consistent with 
the superposition of orthogonal polarization modes (OPMs), and this only at 328 MHz and in superposition with an apparently 
randomly polarized component. The observed periodic pattern is decomposed into the sum of two orthogonally polarized, out- 
of-phase drift patterns, one of which shows a dramatic jump in subpulse phase near the leading edge of the pulse window, which 
probably relates to the phase jump earlier reported in total intensity at 1380 MHz. For PSR B0320+39 and PSR B0818-13, 
considerable periodic fluctuations away from OPM orientations are seen, a condition that also occurs in the trailing half of the 
pulse in PSR B0809+74 at 1380 MHz. In some cases the deviation is so strong that the periodic locus of the polarization vector 
in the Poincare sphere is almost circular, in contrast to the strictly colinear states of superposed OPMs. Several possibilities are 
discussed for the physical origin of these patterns. The similarity between the subpulse patterns in one of the OPMs of PSR 
B0809+74 at 328 MHz to that of the total intensity signal at 1380 MHz supports a picture of superposed, out of phase drift 
patterns. To explain the full range of behaviour seen in the three pulsars, it must be possible to produce at least three arbitrarily 
polarized superposed patterns. While the data do not suggest a particular approach for the empirical decomposition of patterns 
into non-orthogonally polarized components, the specific, quantitative nature of the results should provide strong constraints 
for theoretically driven modelling. 
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1. Introduction 

In the analysis of the bewildering variety of shapes and polar- 
izations of individual pulses, one tends to seek order, in the 
hope that an understanding of simple, repeatable phenomena 
might contribute to a more complete picture of the factors in- 
volved in pulsar emission. In this work two such phenomena 
are considered: the tendency of polarization states to cluster 
around two orthogonal orientations (orthogonal polarization 
modes; OPMs) and the periodic occurrence of subpulses that 
"drift" across the pulse window (drifting subpulses). Both phe- 
nomena offer potential insight into the underlying physics, be- 
cause they can be tied by simple models to the geometry of 
viewing emission that is beamed along dipolar magnetic field 
lines, if the polariz ation is oriented parallel or perpendicular to 
the m a gnetic field jRadhakrishnan & Cookell96^lKomesarofil 



1 19701 Manchester etalJ 1 1975t backer et all 1 19761) . and "the 
drifting subpulses correspond to emitting entities that r otate as 
an ensembl e about the magne tic axis jRudermanl 1 972t IWrightl 
ll98lHAshworthll988tbdwards & Stappersl2002l) . 

That the polarization state is related to the subpulse mod- 
ulation pattern has been known for some time. Due to insuf- 
ficient instrumental time resolution, it was initially thought 
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that a smooth transition of position angle occurred over the 
course of each subpulse dTavlor et al.lll 97 ll) . It was later con- 
clusively shown using data from four pulsars with drifting sub- 
pulses, that the changes are generally sharp jumps from one 
orthogonal state to another, with the point of transition be- 
ing relat ed to the subpulse pos ition rather than the average 
profile (Mancheste r et alJll975l) . Two of these pulsars (PSR 
B0031-07 and PSR B0809+74) show quasi-periodic drifting, 
with the consequence that, for any given pulse longitude, the 
OPM state alternates periodically at the period of the subpulse 
drifting. Using a high-resolution digital pulsar backend and 
modern v isualisation techno logy, this picture was confirmed by 
iRamachandran et al.ld2002h in observations of PSR B0809+74 
at 328 MHz. A natural interpretation of this phenomenon in 
terms of the sta ndard rotating spa rk "carousel" model of drift- 
ing subpulses (Ruderman 1972) is that the observed beam 
is the incoherent superposition of two orthogonally polar- 
ized, azimuthally offset carouse l beams fedwards et alJ l2003: 
iRankin & Ramachandranl200 3). however to date a quantitative 
test of the model and determination of the required azimuthal 
offset have not been performed. 

Much empirical model-building has been done on the 

basis of the interpretation of data in terms of the canon- 

I — 1 1 1 

ical geometric models of polarization (e.g. Rankin 1983; 
iLvne & ManchesteJfl988l lRankinlfl993l iMitra & D eshpande 
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1999; H n & Manchester! l200ll) and drifting subpuls e s (e.g. 
Deshpande & Rankinl 1 19991 beshpande & Rankinl | 200l| 
Asgekar & Deshpandd 1200 it Ivan Leeuwen et alJ ho03[ 
Gil & Sendvkll2003URankin et al.ll2 003). n practice, however, 
significant deviations from the models are frequently observed. 
In many pulsars, the distribution of polarization position 
angles in individual pulses is broadened beyon d what is 
expec ted from the instrumental noise (e.g. Stine bring et alJ 
1984), and a distribution with bimodal peaks separated by an 
angle other than 90° is often observed (e.g. Backer & Rankin 
1980). Related to this, mean polarization p rofiles are generally 
poorly fit by the simple geometric model dEverett & Weisberg 
120011). Empirical models involving non-orthogonal modes 
jMcKinnonll2003l) and superposed randomly polarized radia- 
tion (RPR: lMcKinnonl2 004) have been pursued as a first-order 
explanation for non-OPM be haviour, however in at least one 
known case (PSR B0329+54: [Edwards & StapperslEooH) . the 
distribution of polarization states on the Poincare sphere shows 
a remarkably complicated form indicating a more complex 
origin, most likely in propagation effects. Likewise, in the 
case of drifting subpulses, deviations from the predictions 
of simple geometry appear to be the rule rather than the 
exception dEdwards & Stappersll2003bi) . It is important that 
these deviations are characterised and understood if physical 
information is to be garnered from observations of polarization 
and drifting subpulses. 

In considering the deviations from model prediction s in the 
shape of drift bands in PSR B0320+39. Edwards et all d2003l) 
arrived at the model of offset carousel beams mentioned ear- 
lier. They also noted that, while the beams should be roughly 
orthogonally polarized for PSR B0809+74, to explain the to- 
tal intensity observations of PSR B0320+39 the beams must in 
this case be incompletely and/or non-orthogonally polarized. 
This is due to the fact that the two beams, offset by half the 
spacing of subbeams, must be of equal mean intensity near 
the centre of the pulse window in order to explain the com- 
plete lack of periodic subpulse modulation seen there through 
destructive interference. The superposition of equal amounts 
of completely, orthogonally polarized radiation is incompatible 
with the strong polarization present in this region of the mean 
pulse profile. Therefore, the true origin of polarized drifting 
subpulse patterns must in general be more complicated than, or 
at least different to, the model o f offset OPMs considered by 
iRankin & Ramachandrarjj d2003l) . 

In this work we present an analysis of the polarization fluc- 
tuations associated with periodic subpulse modulation in PSR 
B0320+39, PSR B0818-13 and PSR B0809+74, with a view 
to quantifying not only the OPM behaviour and the associated 
phase offset between the modes, but also deviations from OPM. 
The observations and analysis techniques are described in Sect. 
|2] Sect. [5] presents the observational results, beginning with 
PSR B0809+74 as its behaviour proves easiest to understand. 
This is followed by the discussion (Sect. |4j and conclusions 
(Sect.|5}. 



2. Observations and Analysis 

2.1. Observations 

The observations used in this work were taken from the data 
archi ve accumulated sin ce 1999 with the PuMa pulsar back- 
end dVoute et al.l 120021) of the Westerbork Synthesis Radio 
Telescope. Of those pulsars known to show nearly coher- 
ent periodic drifting subpu lses, observatio ns were available 
of PSR B0320+39 dlzvekova et alJ 1 19821) . PSR B0818-13 
dRitchingsll976tlLvne & AshworfrJ 19831) and PSR B0809+74 
( Sut ton et all 1 970) in the 92 cm wavelength band, and for PSR 
B0809+74 also in the 20 cm band. At the start of each ob- 
serving session, a linearly polarized point source was observed 
with the interferometric correlator, and the relative phases of 
both polarization channels of each of 14 telescopes deter- 
mined. Given these phases, and the location of the source to 
be observed, the signals were added in each polarization of 
1-8 10 MHz-wide observing bands, and the resultant tied ar- 
ray signals processed by PuMa in its digital filterbank mode, 
producing the four Stokes parameters as a function of time 
and frequency. A full polarization calibration is not possi- 
ble due to non-linear system components, however the qual- 
ity of polarization measurements is sufficient if the grossest 
source of potential error, the relative phase of system gain 
in the X and Y polarizati ons, is determined and corrected 
dEdwards & Stappersll2004t) . In principle the phase difference 
is removed as a result of the telescope phase alignment per- 
formed at the start of each observing session, however for 
unknown reasons this often fails for observations in the 92 
cm band, which can be noticed by the presence of Faraday 
modulation in the nominal Stokes V, and by comparison with 
previously published polarimetric pulse profiles. For this rea- 
son, observations in this band were calibrated for differen- 
tia l gain and phase in o ffline processing using the technique 
of [Edwards & Stappers (2004). Samples were then combined 
across all frequency channels, including compensation for the 
effects of interstellar dispersion and Faraday rotation, and the 
resultant time series were divided into segments equal in length 
to the apparent pulse period, to form two dimensional arrays in 
pulse longitude and pulse number. It is important to note, par- 
ticularly in the case of PSR B0809+74 where observations at 
two frequencies were used, that the alignment of "pulse longi- 
tude" to the rotational phase of the pulsar includes an arbitrary 
offset. 

2.2. Analysis 

For the purposes of qualitative visualisation and comparison 
with earlier works, use is made use of plots of individual pulse 
sequences, mean drift bands (formed by folding the data mod- 
ulo the subpulse modulation period), and histograms of the po- 
larization orientation. However, these techniques do not extract 
all of the available information about polarization fluctuations, 
and to this end we define and use techniques that are a general- 
isation of the longitude-resolved fluctuation spectrum (LRFS; 
I — II — ~i 

Backer 1970) to the four Stokes parameters. These techniques 
allow for a more sensitive and quantitative characterisation of 
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the polarization fluctuations. Along with a general description 
of the form of fluctuations expected under orthogonal polar- 
ization modes (OPMs), polarime tric functions a nalogous to the 
LRFS were presented by lEdwards & Stappersl J2004I) for the 
limited purpose of removing the effects of scintillation on ape- 
riodic fluctuation statistics. The descriptions are repeated here 
with a particular view to application of these techniques to 
quasi-periodic polarization fluctuations. 

The (complex) LRFS is obtained by taking the one- 
dimensional discrete Fourier transforms along vectors of con- 
stant pulse longitude: 

where i = y—l, N is the number of pulses, Ip is the total inten- 
sity (Stokes I) in pulse longitude bin / of pulse number j, and 
angle brackets denote averaging over pulse number. The cor- 
responding information for fluctuations in the abstract vector 
space of the Poincare sphere, pji — (Q U V) T , is given by the 
"polarization LRFS" (PLRFS): 

P k l = ^Yje-^ N ( Pjl -( Pl )). (2) 

7=1 

Any sequence of N integer values of k provides complete infor- 
mation on the fluctuation statistics, given the aliasing effected 
by the finite sample interval (one pulse period). It is most con- 
venient to take k in the range (—N/2,N/2], in the knowledge 
that the part of the spectrum with k > corresponds to the 
spectrum of the analytic signal, while the fact that the input 
signal is real-valued gives the result that Pu = P* . 

{—k)l 

As with the complex LRFS, the PLRFS decomposes the 
signal in each longitude bin into a sum of sinusoids. Because 
subpulse modulations are not perfectly periodic, the fundamen- 
tal response is not confined to a single coefficient of the spec- 
trum, and some means to account for this is necessary for op- 
timal sensitivity. The approach taken here is to take the LRFS 
and PLRFS of short (128-pulse) segments of data, over which 
the response is typically confined to one coefficient, and to 
add together, for each longitude bin, the appropriate coefficient 
from each of the spectra, after compensating them for an ar- 
bitrary, longitude-independent phase offset. The latte r is deter- 
mined iteratively using the algorithm given by Edw ards et alJ 
l l2003l) . The results are longitude-dependent complex envelopes 
in total intensity and polarization. Because the phase and am- 
plitudes can be different in each of the three components of 
p, elements of the PLRFS or the complex polarization enve- 
lope in ge neral describe ellipses in p-space (see also, for exam- 
ple, |Bon^^Wolij[l99Sj for a general treatment of complex 3- 
vectors as phased ellipses). The parameters of the ellipse may 
be obtained by observing that any complex vector P may be 
written in terms of the real orthogonal vectors describing the 
semi-major and semi-minor axes A and B and the phase (cp) at 
which P is measured 

P = (A + iB)e i *. (3) 
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The phase may be determined (modulo n) from the argument 
of P 2 = P P: 

P 2 = (A 2 - B 2 ) e 124 , (4) 

and used to obtain ±A and +B as the real and imaginary parts 
of Pe^. Since orthogonal polarization states are anti-parallel 
in /7-space, fluctuations due to pure OPM behaviour should be 
strictly colinear, i.e. \B\ = 0. In this case, the phases and am- 
plitudes of modulations in the two modes can be determined as 
the arguments and moduli of m\ and ma, from the relations 

L = m\ + ni2 (5) 
p = Ae® = l A (mi - m 2 ) , (6) 

where 1 A = A/\A\ is a unit vector parallel to the polarization 
orientation of mode 1 (i.e. the direction in which all fluctua- 
tions occur), L and P are elements of the LRFS and PLRFS or 
of the determined complex envelopes, and the longitude depen- 
dence is dropped for clarity. These relations describe the inco- 
herent superposition of radiation in two orthogonally polarized 
modes, and can be trivially solved to give 

mi = (L + P ■ 1 A ) 12 (7) 
m 2 = (L-P-W/2. (8) 

Note that the PLRFS is the analogue of the complex LRFS. 
In situations where phase information is irrelevant, or where the 
pulse sequence is to be divided into segments and the resulting 
spectra added (e.g. if the desired spectral resolution is much 
lower than the number of pulses available) , LRF power spec- 
tra are generally used (e.g. lBackerll 19731 Edwards & Stappers 
l2003al) . Unlike the LRFS, for the PLRFS the analogue of the 
power spectrum is not simply obtained by taking the element- 
wise squared magnitudes of the complex spectrum, for this 
would neglect possible correlation between the three compo- 
nents of p. The complete autocorrelation statistics of p must 
also include covariance terms, a condition satisfied by the spec- 
tral density tensor 

S kt = PuPu\ (9) 

where ' denotes the Hermitian transpose. Unlike the PLRFS, 
the spectral density tensor contains no absolute phase infor- 
mation, and can thus be summed over multiple segments of 
data without the need to track phase from segment to segment 
(which in any case is only possible for nearly coherent signals). 
Viewed as a two-dimensional array of 3 x 3 matrices, its ele- 
ments are the coherence matrices of p, for a given fluctuation 
frequency (k/N cycles per pulse period) in a given pulse longi- 
tude bin (I). 

Summing the tensor over a range of k provides infor- 
mation on the fluctuation statistics of the data over the cor- 
responding frequency range, with the limiting case that the 
sum over all k gives the real-valued covariance matrix of 
the data. In practice, small values of \k\ from are excluded 
the sum to remove bias due to scintillation-induced fluctua- 
tions, and the mean spectral density tensor (estimated from 
off-pulse longitud es) is subtracted to remove measurement 
noise bias (Edwards & Stappers 2004). The covariance matrix 
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can be subject to eigenvector decomposition to provide in- 
formation on the principal directions of th e fluctuations and 
the shape of the cloud for med in p-space jMcKinnonll2004l 
Edw ards & Stappersll2004l) . These are specified in terms of an 
orthogonal basis composed of the three unit eigenvectors, in 
combination with the three eigenvalues that give the variance 
of p along the corresponding directions. The covariance of p 
between components in any pair of eigenvectors is zero, thus 
the decomposition effectively models the distribution of p as 
an ellipsoidal cloud, with the eigenvectors as its axes. In ad- 
dition to the individual eigenvalu es, use shall be made of the 
polarization entropy dEdwards & Stappers 2004), which com- 
bines the information from the three eigenvalues to give a mea- 
sure of how isotropic the fluctuations are, on scale of (all 
points colinear) to 1 (completely isotropic, spherical distribu- 
tion). Since OPMs define antiparallel orientations in />-space, 
the signature of pure OPM behaviour is perfect linear corre- 
lation, with only one significantly positive eigenvalue (with a 
corresponding vector aligned with the modal orientation) and 
an entropy of zero. 

In principle one may also measure the covariance matrix 
summed over k > 0, or a subset of positive k, to obtain 
the coherence matrix of the analytic signal, or a desired fre- 
quency interval of the analytic signal, corresponding for ex- 
ample to quasi-periodic modulations. The result could be de- 
composed using eigenanalysis into the incoherent sum of or- 
thogonal elliptical fluctuations, or analysed us ing other tech- 
niques of three-dimensional pola rimetry (e.g. ISamsonl|l973t 
ICarozzi et alJbOOCt |Pennisl l2004). However, the pulsars cov- 
ered in this work have sufficiently coherent fluctuations that the 
PLRFS and the real (scintillation-corrected) covariance matrix 
provide all of the required information on the periodic and non- 
periodic fluctuations observed. 

3. Results 

3.1. PSRB0809+74 
3.1.1. 328 MHz 

The observations of PSR B080 9+74 pre sented here are the 
same ones used by iRamachand ran et all J2002t) . who found 
that, at 328 MHz, each subpulse experiences an abrupt tran- 
sition between two elliptical OPMs. Unfortunately, it appears 
that this observation was affected by the the failure of the pre- 
observation phase calibration (Sect l2. II . The resultant X-Y rel- 
ative phase offset (which was measured as -51°) causes a rota- 
tion of p about the Stokes Q axis, resulting in mixing of U and 
V. Although this effect is complicated by the attempted correc- 
tion for Faraday rotation made before summing the frequency 
channels, the net effect was that some of the linear polarization 
masqueraded as circular polarization (and vice-versa). When 
correctly calibrate d, the pattern is qual itatively similar to that 
presented by Ramachan dran et aD J2002I) . except that periodic 
sign flips appear completely abs ent in Stokes VI I (Fig. Q3 , in 
agreement with the findings of Manchest er et alJ 41975). The 
second row of Fig. ^ shows the mean drift band, produced 
by folding a sequence of ~ 300 pulses modulo the drift pe- 
riod, and plotting the result in terms of total and polarized in- 



tensity, and polarization orientation (position angle and ellip- 
ticity angle 1 ). This sequence was selected as being unaffected 
by nulls and the associated modulation period variations JColel 
ll970HLvne & Ashworfhll983l:lvan Leeuwen et alJ2003h which 
would smear out the drift band shape. Clearly seen is a reduc- 
tion in the polarized intensity in the middle of the drift band, 
accompanied by sharp jumps in position angle. No analogous 
jump is seen in ellipticity angle, as would be expected for ellip- 
tical OPMs. The longitude-resolved histograms of position and 
ellipticity angles (Fig. |3 confirm this picture, and show that 
major random deviations from the nominal linear OPM orien- 
tations occur in both parameters. 

The results of eigenanalysis of the covariance matrix (Fig. 
provide quantitative support for the features noted above. 
The very small values of ellipticity angle for the dominant 
eigenmode concur with the lack of sign flips in Stokes V// (Fig. 

indicating that the OPMs are very nearly purely linearly 
polarized. Incidentally, we note that such a condition would 
be most unlikely under an arbitrary X-Y relative phase error 
(especially given the swing in position angle with pulse lon- 
gitude), confirming the quality of the calibration. The position 
angle swing of the dominant eigenvector is nearly perfectly lin- 
ear, being much more robust to slight deviations from orthog- 
onality than the position an gle of the average profile (see also 
Edwards & Stappers 2003b). The second and third eigenvalues 
are both significant and nearly equal in value at most pulse lon- 
gitudes, indicating that the distribution in the Poincare sphere 
is roughly that of a prolate spheroid. This is consistent with the 
superposition of an apparentl y randomly polar ized component 
on top of the OPM radiation ( Mc Kinnonl2004 . 

To examine the relationship of OPMs with the drifting sub- 
pulse modulation in more detail, the longitude-dependent am- 
plitude and phase of the quasi-periodic total intensity and po- 
larimetric fluctuations were computed using the methods of 
Sect. 12.21 The result is shown in Fig.|4] The measured values of 
|B| (semi-minor axis) are small, showing that the periodic fluc- 
tuations in p are nearly confined to a line, consistent with an 
origin in OPMs. Working from this assumption, the observed 
periodicity was decomposed into the sum of two out-of-phase 
OPMs using Eqs. [7]and[8] The resultant longitude-dependent 
amplitude and phase envelopes for the two modes are also de- 
picted in Fig.0] along with their phase difference (for conve- 
nience). The amplitude peaks are offset in pulse longitude, and 
for most of the pulse window the phases are offset by ~ 50°. 
However, at pulse longitude ~ 47° there is a sharp "jump" of 
~ 120° in the phase envelope of one of the modes. 

3.1.2. 1380 MHz 

Owing to the reduced depth of subpulse modulation of PSR 
B0809+74 at high frequencies JNowakowski et aljfl 982). it is 
difficult to obtain sensitive measurements of its subpulse prop- 
erties at 1380 MHz. Although the best polarimetric observation 
available was of significantly lower quality than the total inten- 

1 The ellipticity angle ix = \ tan -1 V/(Q 2 + U 2 ) il2 ) complements the 
position angle to completely specify the orientation of a polarization 
vector in QUV space, independent of its length. 
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Fig. 1. Colour plots of polarimetric subpulse patterns, pulse number (ordinat e) versus pulse long itude (abscissa). The top row 
shows a sequence of single pulses PSR B0809+74 at 328 MHz, after Ramachandran et al. (2002) but now correctly calibrated. 
The other rows show mean drift bands, produced by folding the Stokes parameters modulo subpulse modulation period, and 
plotted twice for continuity. From top to bottom: PSR B0809+74 at 328 MHz, PSR B0320+39 at 357 MHz and PSR B0818-13 
at 328 MHz. In the folded plots, a parameterisation of the polarization that splits it into intensity (\p\) and orientation (position 
angle, and ellipticity angle, x - \ tan -1 V/(Q 2 + U 2 ) 1 ^ 2 : see footnote to Sect. B.l.H are used to better assess consistency with 
OPMs. The colour scales used are shown at the bottom. Intensity parameters (left colour scale) are normalised by the peak value 
plotted in a given panel. 
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Fig. 2. Polarization histograms. The top panel shows the total intensity (solid) and polarized intensity (dotted) of the average pulse 
profile. The centre panel shows the longitude-dependent histogram of position angle (plotted twice for continuity) along with the 
position angle of the mean profile. The bottom panel shows the histogram of sin 2% (where \ is the ellipticity angle), along with 
its counterpart from the mean profile, (sin 2^- is chosen over^- so that each bin samples an equal area of the Poincare sphere; 
see Edwards & Stappersl l2004.) For convenience the scale is drawn (non-linearly) in terms x- Only samples with a polarized 
intensity greater than three times the RMS noise (estimated from the off-pulse) were included. The density scale is such that bins 
containing no samples are white, while other values are mapped linearly from medium grey to black. 



sity observation used by Edwards & Stappersl J2003bl) (which 
apparently encountered extraordinarily strong amplification by 
the interstellar medium), some information on the subpulse be- 
haviour can nevertheless be obtained. Th e position angle his- 
togram is similar to the scatter plot of iRamachandran et alJ 
(2002), while the ellipticity angle histogram is broad and fea- 
tureless, so for brevity these are not displayed. 

Turning to the results of eigenanalysis (Fig. [3}, we see that 
the behaviour is very different on either side of pulse longitude 
~ 38°. On the leading side, only one significant eigenvalue is 
measured, indicating very pure OPM-associated fluctuations. 
On the trailing side the fluctuations become almost isotropic, 
accompanied also by a ~ 90° step in the polarization position 
angle of the mean profile. 

Further insight is given by the decomposition of the peri- 
odic polarimetric fluctuations into their semi-major and semi- 
minor axes (Fig. |4}. Firstly, it is seen that on the leading part 
of the profile, the periodic fluctuations are approximately lin- 
ear (and of similar amplitude to those in Stokes I), while in the 
trailing part significant values of \B\ are detected, indicating 



that the pattern cannot be produced by the superposition of two 
OPMs. The phase envelope of the semi-major axis (A) closely 
follows that of the total intensity in the leading part of the pro- 
file, consistent (along with the equal amplitudes) with the pres- 
ence of modulations in only one of the polarization modes. The 
step of ~ 120° in th e phase envelope of th e total intensity (pre- 
viously reported by Edwards & Stappers 2003b) is accompa- 
nied by a smoother, smaller transition of the opposite sense in 
the phase of the semi-major axis, although the meaning of this 
is not clear in the absence of a model for the origin of elliptical 
periodicities. 

3.2. PSR B0320+39 

The orientation histograms (Fig. for PSR B0320+39 at 
357 MHz show broad distributions in position angle and es- 
pecially in ellipticity angle. This character is confirmed by 
the results of eigenanalysis (Fig.|5|l, which show that all three 
eigenvalues are generally significant, and all different in value. 
The position angle of the primary direction of fluctuation dif- 
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PSR B0809+74 (328 MHz) PSR B0809+74 ( 1 380 MHz) PSR B0320+39 (357 MHz) 



PSRB0818-13 (328 MHz) 
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Fig. 3. Results of eigenanalysis of the covariance matrices. The top panel shows the average intensity profile (thick solid line), 
the square roots of the eigenvalues (solid, dashed, dotted thin lines, in descending order of value), and the polarization entropy 
(thick dotted line). The middle panel shows the position angle of the mean polarization vector (thick line, plotted repeatedly at 
offsets of 180° for continuity) and of the first eigenvector (thin line, plotted repeatedly at offsets of 90° for continuity under the 
equivalence of antiparallel vectors as eigenvectors). The bottom panel shows the ellipticity angle of the mean polarization vector 
(thick line) along with the ellipticity angle of the first eigenvector, plotted twice with opposite signs. 



fers markedly from that of the average profile. The amplitude 
and phase envelopes of the periodic polarization fluctuations 
(Fig. 0J| reveal, via the significant values for \B\, that fluctu- 
ations in at least two directions in QUV space contribute to 
the observed periodic patterns. In fact, this can be seen in the 
mean drift band (Fig.Q: around pulse longitudes 17-18° and 
20-21° the periodicity in ellipticity angle is clearly smoothly 
varying, in contrast with the sign-reversing rectangular wave- 
form expected under OPM. There is some evidence for similar 
behaviour in the position angle around longitude 18-20°. Fig. 
shows that semi-major axis fluctuates roughly in phase with 
the total intensity, including the ~ 180° step in phase (longi- 
tude ~ 18.5°), which more sensitive observations have shown 
is very sharp and accompanied by an al most complete lack o f 
periodic modulations in total intensity (Edwar ds et alj E003). 
In contrast to this, significant periodic polarization fluctuations 
(about an elliptical locus) are detected throughout this part of 
the pulse window. 



derstanding the strange behaviour around longitude 39.5° is 
given in the distribution of ellipticity angles : most points scat- 
ter about the V/\p\ = 1 pole of the Poincare sphere. This is con- 
firmed by Fig. [51 which shows the two-dimensional distribution 
on the Poincare sphere in this longitude region. Although the 
position angle distribution is roughly bimodal, the true two- 
dimensional distribution is unimodal, with an extended shape 
that prefers some position angles over others. The results of 
eigenanalysis (Fig.[3ji confirm the chaotic picture given by the 
histograms, which taken together show that there is no partic- 
ular reason to consider OPMs to play a role in the polarization 
fluctuations. The periodic part of the fluctuations (Fig. |4jl are 
just as complicated, with the pattern tracing an elliptical locus 
in p-space, which is almost circular at longitude ~ 39.5° (as 
indicated by the near equality of \A\ and \B). There is also ev- 
idence for this in the mean drift band (Fig. 0, in the form of 
continuous variations in the ellipticity angle, and roughly bi- 
modal variations in the position angle. 



3.3. PSRB0818-13 

The position angle histograms of PSR B0818-13 at 328 MHz 
(Fig. [3 reveal a chaotic jumble of preferred polarization direc- 
tions as a function of pulse longitude. The distributions clearly 
differ strongly from the expectations of OPM. A clue to un- 



4. Discussion 

The observations presented here prove what was already sus- 
pected in the introduction: that periodic polarization modula- 
tions are more complicated than anything that can be produced 
by the out of phase superposition of two orthogonally polarized 




Fig. 4. Amplitude (top) and phase (bottom) of subpulse modulations, in total intensity (thick solid line) and polarization (semi- 
major axis (\A\): thick dashed, semi-minor axis QUI): thick dotted). The phase of the semi-minor axis is not plotted, since by 
definition it is offset from that of the semi-major axis by ±90° In the left-most column only, the amplitude and phase of two 
OPMs (thin dotted and dashed lines) as inferred from the total intensity and semi-major axis values are plotted along with the 
phase difference (thin solid line). A constant slope was subtracted from all phases for clarity. The (dimensionless) slopes used 
were 25 for PSR B0809+74, -60 for PSR B0320+39, and 60 for PSR B0818-13. 




Fig. 5. Distribution of polarization orientations in pulse longi- 
tudes 39-40° for PSR B0818-13 at 328 MHz. The projection 
is Lambert's azimuthal equal area, interrupted at the equator 
and centr ed on the poles of the Poincare sphere (V/|p| = +1) : 
see lEdwards & S tappers! J2004I) for details. Meridians and par- 
allels are marked at intervals of 30°, with the upward vertical 
lines corresponding to a position angle of 0°. 



drift patterns (cf. lRankin & Ramachandranl2003l) . Of the three 
pulsars studied, only PSR B0809+74 fits this picture, and this 
only at 328 MHz. This is subject to the caveat that a significant 
component of randomly polarized radiation (RPR) appears to 
be present. 

It seems pertinent to ask whether some other mecha- 
nism is responsible for the periodic polarization fluctuations 
seen in pulsars with drifting subpulses. Whatever the under- 
lying mechanism for the production of drifting subpulses, it 
is somehow capable of modulating the total intensity of the 



emission produced in the magnetosphere. The local magne- 
tospheric conditions also determine the path taken by plasma 
wave "rays" in their propagation through the magnetosphere 
(e.g. Russell & Luo 2004; IPetrovalEoOOl) . the alterations im- 
posed on their polarization through birefringent effects in the 
polarization limiting region (e.g. ICheng & Rudermanl Il979t 
iLvubarskii & Petroval IT999). and, potentially, the d egree of 
conve rsion to an orthogonal sense of polarization JPetroval 
1200 ill . It therefore seems plausible that these processes could 
also be modulated at the subpulse fluctuation period, giving rise 
to periodic changes in the polarization state. Only through care- 
ful modelling could the consistency of the observations with 
these effects be checked. 

However, there is a striking feature of the 328 MHz ob- 
servations that is particularly suggestive that the modulations 
do indeed arise due to the out of phase superposition of two 
or more periodic patterns. When the measured fluctuations 
were decomposed as the sum of out of phase OPM mod- 
ulations, the longitude-dependent phase of one of the two 
modes was seen to execute a sudden jump of 120°, accompa- 
nied by a reduction in the amplitude of fluctuations. A very 
similar featur e is seen in the total inte nsity modulations at 
1380 MHz llEdward s & Stappers 2003b), which, as shown in 
this study, occur in a single polarization mode (at least in the 
leading part of the pulse window). In our view, such a corre- 
spondence is highly unlikely to arise unless the decomposi- 
tion performed at 328 MHz has a true phy sical basis. As al- 
ready argued by Edwards & Stappers (2003b), the behaviour at 
1380 MHz could be explained via superposed subpulse patterns 
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with longitude-dependent phase offsets of an arbitrary amount 
(i.e. not just 180°). If the same effect causes a phase jump in a 
single OPM (as observed at 328 MHz), which itself is appar- 
ently accompanied by an out-of-phase subpulse pattern in the 
complementary polarization mode, then there are at least three 
offset subpulse modulation components present in B0809+74. 

A second important implication of the phase jump in one 
OPM of PSR B0809+74 at 328 MHz concerns the sam- 
pling of the polar cap effected in observations at low ver- 
sus high frequencies. The apparent association of this fea- 
ture with the phase jump in total intensity at 1380 MHz 
supports the frequency-dependent pulse profile alignme nt de - 
rived from time-stamped observations by Kuzm inet alJ {l998), 
which differs from the (arb itrary, visual) alignment used by 
lEdwards & Stappersl J2003bl) . Under the updated alignment, 
the leading and trailing edges of the 1380 MHz mean profile 
occur significantly earlier than the corresponding features of 
328 MHz pulse, arguing that observations at either frequency 
samples a region of magnetic azimuth not sampled by the other, 
and/or that symmetry-breaking effects such as abberation, re- 
tardation, refraction or asymmetric or patchy polar cap patterns 
play a role. This should be contrasted to the picture of low- 
frequency "absorption" of emission solely on the leading edge 
of the profile, arrived at bv lBartel et alJ lll98lh . Their result de- 
rived from an unphysical frequency-dependent longitude scal- 
ing performed to expl ain a difference in drift rate th at is actually 
due to the phase step (Edwards & Stappers 20033). 

The polarization fluctuations seen in PSR B0320+39 and 
PSR B0818-13 have a different character, in that the mean lo- 
cus followed by the tip of p under the periodic modulation is 
elliptical, rather than linear as would be expected under the su- 
perposition of OPMs 2 . Therefore, if the picture of superposed, 
out of phase drift patterns is also to explain the modulation in 
these pulsars, the component patterns must be allowed to have 
polarization states that are not orthogonal (or aligned). That 
this should be the case is not surprising, for whatever effect 
is responsible for the offset in modulation phase could plausi- 
bly also affect their polarization, since both are believed to be 
tied to the differences in effective viewing geometry. Indeed, 
such considerations could potentially explain the slight non- 
orthogonality observed even at 328 MHz in PSR B0809+74. 

Further insight into the phenomenon of superposed 
drift modes is likely to come only with a larger sam- 
ple size. Two previous mea surements are of relavance here. 
iDeshpande & Rankinl {2001) found that subbeam patterns in 
two OPM-segregated polar cap images of PSR B0943+10 
were offset by roughly half the subbeam spacing. Similarly, 
iRankin & Ramachandranl (|2003) found that the phases of the 
fundamentals in the longitude-resolved fluctuation spectra of 
the two OPMs of PSR B 1237+25 were offset by roughly 180°. 
Unfortunately, the disjoint mode segregation algorithm used in 
both measurements is prone to produce anti-phase modulation 



2 Purely sinusoidal modulations are confined to elliptical paths by 
definition, but the presence of higher harmonics, neglected in our 
analysis, could in principle allow for any closed locus in p space. 
Nevertheless, the ellipticity of the fundamental is proof that the locus 
is not confined to a line. 



(d) 




Fig. 6. Illustration of the results of three-w ay polarization mode 
segregation (Deshpande & Rankin 2001) on the superposition 
of sinusoidal OPMs. (a) Simulated intensity of two OPMs (in a 
given longitude bin) as a function of time, (b) Resultant signal 
in the component of p corresponding to the modal orientations, 
(c) Primary polarization mode intensity signal as inferred by 
three-way segregation (solid) and its fundamental Fourier com- 
ponent (dashed), (d) As for (c) but for secondary mode. Note 
that the underlying OPMs are nearly in phase, while the funda- 
mentals of the inferred signals are in pure anti-phase. 



even if the patterns arise from the nearly in-phase superposi- 
tion of OPMs, since the modal patterns cannot overlap in the 
segregated signal. It is easily shown (Fig.|6} that in the patho- 
logical case of superposed sinusoidal modulations, the period- 
icity will be in complete anti-phase between the two inferred 
modes, regardless of the actual phase offset. For an accurate 
measurement, firstly it must ascertained whether the data are 
indeed quanitatively consistent with pure OPM, and secondly 
the mode separation method must take the superposition of ra- 
diation into account. 

Unrelated to the phenomenon of drifting subpulses, the 
1380 MHz results for PSR B0809+74 allow for some inter- 
esting constraints to be placed on the origin of randomly po- 
larized radiation (RPR), in light of the fact that it only oc- 
curs in the trailing half of the pulse window (Fig. [3}- That this 
should b e accompanied by sudden, almost complete d epolar- 
ization JMorris et alJll98lt Ikamachan dran et alJl2002h can be 
understood in one of two ways. One possibility is that the de- 
polarization occurs due to the superposition of two OPMs of 
nearly equal intensity, with the RPR being associated with only 
one of the modes. The second possibility is that the RPR is not 
a superposed emission component, but rather arises due to a 
propagation effect, randomising the polarization of incoming 
radiation and causing depolarization of the average profile. It 
should be noted that this is different to the types of depolaris- 
ing propagation effects usually encountered, where the electric 
field vector is altered on timescales comparable to the inverse 
of the observing frequency. 

The signature of RPR is significant polarization at an ori- 
entation that is consistent over some time scale longer than the 
averaging interval of the individual data samples, yet randomly 
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varying on some longer time scale (in our case, the pulse pe- 
riod). In the case that the RPR is due to a propagation effect, 
this would imply a variable transfer function (Jones or Mueller 
matrix) with a decorre lationtimescale tha t lies between the two 
timescales mentioned. McKinnon (2004) also mentions the vi- 
ability of RPR as a depolarization mechanism, however he ap- 
pears to assume that it must arise as superposed emission. Such 
a condition is not capable of reducing the polarized intensity, 
as is required for PSR B0809+74. Rather, just like unpolarized 
radiation, it can only reduce the fractional polarization, by in- 
creasing the total intensi ty without add ing polarized intensity to 
the average profile. [McKinnoiJ J2004) considers the possibility 
of stochastic Faraday rotation in the pulsar magnetosphere as 
an origin of RPR, but rejects it because the magnetic field is 
unlikely to vary on short time scales, and in any case the effect 
cannot alter the Stokes V component of p. However, this is an 
inadequate basis on which to rule out an origin for RPR in prop- 
agation effects: birefringent effects in the magnetosphere can 
alter all of the Stokes parameters, as the propagation modes are 
believed to be linearly or elliptically polarized (rather than cir- 
cularly, as is the case for Faraday rotation), and the fluctuations 
could be due to variations in plas ma density and/ or propagation 
path, rather than magnetic field (Petrova 2001 and references 
therein). 

5. Summary and Conclusions 

We have presented a detailed study of the average periodic po- 
larization fluctuations in three pulsars with drifting subpulses. 
The basic result of periodic switching between two OPMs in 
PSR B0809+74 iMancheste r et alJI 197.4 iRamachandran et alJ 
120021) is confirmed at 328 MHz, albeit in superposition with an 
apparently randomly polarized component. The periodic pat- 
tern can be understood as the superposition of drift patterns 
in two orthogonal polarization modes that are offset by ~ 50° 
in subpulse phase. The drift bands are roughly linear, how- 
ever for one of the modes a sudden jump of 120°, accompa- 
nied by a reduction in the amplitude of modulations, is seen 
near the leading edge of the pulse profile. Very similar be- 
haviour is seen in the total in tensity modulations at 1380 MHz 
llEdwards & Stappersl2003bh . which are probably composed of 
a single polarization mode. That such a correspondence is ob- 
served strongly supports the physical validity of the decom- 
position of the pattern as the superposition of out of phase 
drift patterns. Aligning the pulse profiles on the basis of the 
phase jum p agrees w ith results derived from time-stamped ob- 
servations ( K uzmin et alJI 19981) . implying either that different 
regions of magnetic azimuth are sampled at the two frequen- 
cies, with neither frequency presenting a complete image of the 
polar cap, or that other effects such as a non-dipolar magnetic 
field, refraction, aberration and/or retardation are important. 

The results from PSR B0320+39 and PSR B0818-13 re- 
veal that the behaviour deviates considerably from the predic- 
tions of OPM, tracing in general a quasi-elliptical path in the 
Poincare sphere with each cycle of the subpulse modulation. In 
PSR B0818-13, the unimodal, asymmetric clustering of states 
around the V/\p\ — 1 pole of the Poincare sphere leads to a 
deceptively complicated position angle histogram. 



In order to explain the patterns seen in the three pulsars 
studied, considerable complications to the standard models are 
required, with in general no less than three offset, arbitrar- 
ily polarized drift patterns seen in superposition. For the non- 
orthogonal oscillations, no particular empirical basis for de- 
composition is suggested by the data, however it is hoped that 
the specific, quantitative nature of the results will prove highly 
constraining to theoretically-driven models. 
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